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Magnetic monopoles from gauge theory phase transitions 
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Thermal fluctuations of the gauge field lead to monopole formation at the grand unified phase 
transition in the early universe, even if the transition is merely a smooth crossover. The dependence 
of the produced monopole density on various parameters is qualitatively different from theories with 
global symmetries, and the monopoles have a positive correlation at short distances. The number 
density of monopoles may be suppressed if the grand unified symmetry is only restored for a short 
time by, for instance, non-thermal symmetry restoration after preheating. 

PACS numbers: PACS: 11.15.Ex, 11.27.+d, 74.60.Gc 



It is a generic property of grand unified theories 
(GUTs) that magnetic monopoles of mass of the order 
tom ~ 10 16 GeV exist [HQ, and these monopoles would 
have been produced in large numbers in the GUT phase 
transition at Tqut ~ niM 0- Afterwards, pair annihi- 
lations can decrease the monopole density, but estimates 
show that the number density would still be compara- 
ble to baryons [j]. Because the monopoles are 10 16 times 
heavier than protons, this would have caused the universe 
to collapse under its own weight long ago. 

This monopole problem, alongside with several other 
cosmological puzzles, was wiped away by the theory of 
inflation [Bj], as the monopole density would have been 
diluted to a negligible level by a period of accelerating 
expansion. For this to solve the problem, the reheat tem- 
perature at which the universe thermalizes must be lower 
than Tqut- These constraints are even stronger in mod- 
els with non-perturbative effects such as preheating 
since the GUT symmetry can be temporarily restored 
0, and topological defects formed even if the reheat 
temperature is well below Tqut |9tllC(. It is therefore im- 
portant to understand how monopoles are formed to es- 
timate how strong the bounds imposed by the monopole 
problem really are. 

In this paper, I will discuss monopole formation at a 
phase transition that starts from complete thermal equi- 
librium. It is clear that this is not actually the case for the 
GUT transition, because of the high expansion rate and 
the non-equilibrium effects mentioned above. Neverthe- 
less, the assumption of thermal equilibrium simplifies the 
problem significantly and makes it possible to identify the 
physical mechanisms that are responsible for monopole 
formation. Once these mechanisms are understood, their 
effects can be studied in more realistic non-equilibrium 
settings. 

The symmetry broken at the GUT phase transition is 
a local gauge invariance, whereas most of the existing 
literature on monopole formation implicitly assumes a 
breakdown of a global symmetry. The Kibble (or Kibble- 
Zurek) mechanism 0, 1 1 1| , which forms the monopoles in 
the global case, is ultimately based on the observation 
that the direction of the order parameter cannot be cor- 
related at infinitely long distances. Because the direction 
of the order parameter is not gauge invariant, this argu- 



ment cannot be used in GUTs. 

Moreover, gauge symmetries cannot be spontaneously 
broken |12| . For rather generic parameter values, there is 
actually no phase transition at all, but simply a smooth 
crossover between the phases 0|. Does this mean 
that the whole evolution could be adiabatic and thereby 
there would be no monopole formation at all? 

I will present an argument that shows that monopoles 
are still formed. This result is based on causality and 
the conservation of magnetic charge. In fact, Weinberg 
and Lee [TEL Il6| have used somewhat similar reasoning 
to constrain later annihilations of monopoles after the 
phase transition in the context of the Kibble mechanism. 
As I will argue, there are long-wavelength thermal fluc- 
tuations of the magnetic charge in the symmetric phase, 
and they will freeze out forming monopoles. These fluc- 
tuations are physical and well defined, because the high 
and low temperature phases are smoothly connected. As 
there is, in this sense, a high density of monopoles above 
the transition, one could say that we are describing anni- 
hilation rather than formation of monopoles. This is ob- 
viously a matter of taste, but in any case, the monopoles 
do not correspond to localized energy concentrations in 
the symmetric phase and cannot therefore be thought of 
as particles. 

The mechanism presented in this paper is physically 
different from global theories. Both involve a freeze-out 
of long-wavelength degrees of freedom, but in global the- 
ories this happens when the scalar correlation length di- 
verges at the transition point. In our case, everything is 
finite at the (approximate) transition point but the mag- 
netic screening length diverges in the zero-temperature 
limit. The reason why this freeze-out leads to monopole 
formation is also different. As we shall see, the monopoles 
formed in a gauge theory have positive correlations at 
short distances, which is the opposite of what the Kib- 
ble mechanism predicts ITtI . The number density of 
monopoles will also be qualitatively different from the 
Kibble mechanism. 

Let us start by briefly reviewing the standard Kibble 
mechanism 3] , which is valid in the case of global symme- 
tries. For simplicity, we shall discuss the SU(2) symme- 
try group only, but the same arguments should apply to 
SU(5), SO(10) or other possible GUTs. The Lagrangian 



2 



of the global SU(2) scalar theory is 



D„ 



C = Trd^d^ 



2 Tr$ 2 -A(Tr$ 2 ) : 



(1) 



C 



-TrF^F^ + Tr[D^, $] [D», $] 



where $ is in the adjoint representation, and we are as- 
suming that at zero temperature, the SU(2) symmetry is 
broken. To leading order, this means m 2 < 0. 

We shall consider this theory at a non-zero tempera- 
ture T. When the temperature is high enough, the SU(2) 
symmetry is unbroken. We are asking what happens if 
we start from thermal equilibrium in the symmetric phase 
and gradually decrease the temperature so that the sym- 
metry gets broken. As long as couplings are weak, we 
can approximate the equilibrium and near-equilibrium 
dynamics reasonably well by a classical field theor y w ith a 
temperature-dependent mass term m 2 (T) [lMll^l2fll21j . 
Although it is difficult to make this approach quantita- 
tively very accurate, it gives us a way of thinking about 
the dynamics in terms of hot classical fields without 
quantum mechanical complications. 

It only takes a small change in temperature near the 
critical temperature T c to cause the phase transition, and 
this effect is mainly due to the effective mass parameter's 
changing from positive to negative. Therefore, we shall 
simply consider keeping T fixed and varying m 2 . 

Let us now discuss the dynamics of the global theory 
JTJ. In the high-temperature phase, the field $ vanishes 
on the average. In the broken phase, it would ideally 
have a non-zero constant value = 4>o, where Tr<&Q = 
(f> 2 /2 = — m 2 /2A > 0, but this would require ordering of 
the field at infinite distances, which cannot be achieved in 
finite time. Instead, the scalar correlation length £ grows 
as the transition point is approached, but freezes out to 
some finite value £, which is determined by the critical 
dynamics of the system |l l| and ultimately limited by 
causality. 

After the transition, we can imagine that the system 
consists of domains of radius £, between which $ is totally 
uncorrelated. At each point where four of these domains 
meet, there is a fixed, non-zero probability that the field 
cannot smoothly interpolate between the domains with- 
out vanishing at a point. This point is a monopole, and 
therefore this scenario predicts a monopole density 



Kibble 



r 



(2) 



Furthermore, there is a strong negative correlation be- 
tween monopoles at short distances ^3 H3 : Imagine a 
sphere centred at a monopole. If the radius is greater 
than £, each point at the sphere is uncorrelated with 
its centre and therefore insensitive to whether there is a 
monopole inside or not. Consequently, the average wind- 
ing number must be zero and there must be an anti- 
monopolc within distance w £ from each monopole. 

Having reviewed the Kibble- Zurek scenario, let us now 
turn our attention to the gauge theory. The gradients 
in the Lagrangian are replaced by covariant derivatives 



-m 2 Tr$ 2 



A (Tr$ 



(3) 



where = (ig) 1 [D fl ,D v }. 

At a high temperature T and at weak coupling g, the 
phase structure of this theory is given to a good ap- 
proximation by a three-dimensional effective theory (20j , 
which depends on two parameters, g (T/T c — I) and the 
ratio of the coupling constants X/g 2 . There is a line of 
first-order transitions at small X/g 2 which ends at 
a second-order point at around X/g 2 « 0.3 |2^|. To sim- 
plify the estimates, we actually assume that X/g 2 » 1, 
which is above the critical value so that the two phases 
are smoothly connected 0] and all correlation lengths 
are finite. In particular, the scalar correlation length is 
microscopic |23) . of order l/g 2 T, and therefore we ig- 
nore scalar fluctuations. Even then, it is still possible to 
find an approximate crossover point, which separates a 
"symmetric" and a "broken" phase, and we will use this 
terminology although it is not quite precise. 

In the broken phase, there is still an unbroken U(l) 
subgroup, and the corresponding magnetic field is given 
by the 't Hooft operator 



n - 1 



Tr$F jk + ^ r -Tr$(D j &)(D k &) 
2ig 



(4) 



where <J? = $y / 2/Tr$ 2 is well-defined almost everywhere. 
In continuum, Bi is sourceless apart from points where $ 
vanishes, and has sources of integer multiples of 47r/g at 
those points. Being sources of the magnetic field, these 
points are quantized magnetic charges, and therefore we 
call them magnetic monopoles, whether or not they re- 
semble the 't Hooft-Polyakov monopole solution P, Q|. 
It is straightforward to see that magnetic charge defined 
in this way is conserved, i.e., the world lines of these 
monopoles cannot end. 

The same construction has also been done on lat- 
tice 0, , and the conservation law of magnetic charge 
is preserved. Even though we will not discuss numerical 
simulations in this paper, this is very important, because 
classical field theory cannot be in thermal equilibrium in 
continuum. Because the monopoles are well-defined and 
stable on lattice, we can consistently talk about them 
and the magnetic field at a non-zero temperature. 

Deep in the broken phase, where vhm > T, we can treat 
the monopoles as point-like particles. Therefore, we have 
the standard expression for the equilibrium monopole 
density 



eq 



(m M T) 3 / 2 exp (- 



T )' 



(5) 



The monopole mass m,M is roughly ~ <f>/g ~ 

(— m 2 / Xg 2 ) 1 / 2 . When m 2 decreases further, tom grows 
rapidly, which suppresses n e ^. 
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If the monopoles did not have long-range interac- 
tions, they would be essentially uncorrected and behave 
very much like the magnetic field in the Abelian Higgs 
model |2(jL |2?| . albeit in three rather than two dimen- 
sions. There is, however, a magnetic Coulomb interac- 
tion between the monopoles, and we shall see that it 
suppresses their production. This interaction gives rise 
to correlations, which are reflected in the screening of the 
magnetic field by the monopoles [Hi. l2c| , in analogy with 
the Debye screening of the electric field. 

We define the magnetic screening length £ B as the de- 
cay rate of the correlator of 't Hooft field strength oper- 
ators Bi. In equilibrium, it is approximately 



£,b = l/m B 




~~i (6) 



where qm = 47r/ ' g is the magnetic charge of a monopolc. 
Correspondingly, if we define 

Pm = V • B, (7) 
the magnetic charge-charge correlator is 

2 

m B - mB \x-^\ 



(Pm(x)pm(v)) ~ q M n M \ S(x-y) 



4-7r|x— y\ 



Using Eq. J5J, we find that the equilibrium screening 
length behaves as 



( 9 4 



\ Tvn? 



M 



1/4 



i M /2T 



(9) 



Because there is no phase transition, the correlators of Bi 
or pm cannot change qualitatively when we move to the 
symmetric phase. Otherwise, their behaviour could be 
used to distinguish between the phases. Consequently, 
the screening length is always well defined, and we can 
actually use Eq. © to define the monopole density um in 
the symmetric phase. Furthermore, we expect that above 
the crossover, £b ~ {g 2 T)~ 1 , because the only relevant 
scale for equal-time correlations is g 2 T [l^. 1 

If m 2 is decreased at a constant rate, £b would have 
to grow exponentially fast to stay in equilibrium, but, 
obviously, it cannot grow faster than the speed of light. 
In practice, it would grow much slower than this. This 
means that sooner or later the growth rate d^s / dt needed 
for the system to stay in equilibrium exceeds the maxi- 
mum value, and the system falls out of equilibrium. We 
shall denote the time when this happens by t. 



1 The scale g 2 T is generic in hot gauge theories and is known as 
the magnetic screening scale. This magnetic screening refers to 
non-Abelian magnetic fields, and it is uncertain whether it can be 
thought of originating from some kind of monopoles. In contrast, 
the field Bi is, by definition, screened by monopoles, because the 
monopoles were denned as sources of Bi. 



The screening length £b can still keep on growing, but 
so slowly that we can ignore it if we are only interested 
in finding an order-of-magnitude estimate for the initial 
monopole density. Therefore, we define the freeze-out 
screening length £b as £b at the time when it falls out of 
equilibrium, 

£b = &(**). (io) 

At the time of the freeze-out, the monopole density is 

T g 2 T_ 

>B 



n M 



(11) 



Note that by cooling the system slower, we can make £b 
arbitrarily large. When £b ^> (5 2 T) _1 , the typical dis- 



tance d 



-1/3 



between monopoles and antimonopolcs 



is much shorter than the screening length. 

Even after the freeze-out, the monopole density will 
keep on decreasing, but this is now due to pair anni- 
hilations at length scales shorter than £b- These an- 
nihilations smoothen the distribution of monopoles at 
short distances, but they cannot remove them com- 
pletely [3 EfJ. To see this, consider a sphere of ra- 
dius £b ■ The annihilations may reduce the number 
of monopoles inside the sphere to the minimum, but 
they cannot change its net magnetic charge significantly. 
While the net magnetic charge is zero on the average, it 
fluctuates with a root-mean-squared value of 



Qm((b) 



\ 




d 3 xp M {x) 




T£b- (12) 



Since the annihilations cannot reduce the charge below 
this, the monopole density cannot fall below 



_ Qm{£b) _ -l 

n M ~ ~ q M 

Qm^b 




(13) 



We have not shown how to estimate £b, but nevertheless, 
this expression is clearly different from the Kibble-Zurek 
result (j2J) , because of the explicit appearance of g and T. 

Moreover, as long as Qm(£b) ^> qM, there will be 
clusters of monopoles of equal sign, and the number of 
monopoles in each of them can be large if T 3> 
This means that there is a positive correlation between 
monopoles at short distances, very much in the same way 
as in the case of vortices in the Abelian Higgs model |2o . 
|27| and in stark contrast with the Kibble mechanism. 

We can reach the same conclusions by studying the 
time evolution of the magnetic charge correlator in the 
Fourier space. We define the equal-time correlator G(k) 
by 



(pM(k)p M (q)) = q M G(k)(2Tr) 3 6 [k + q 



(14) 
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and from Eq. J5J), we find 

G{k) 



T m 2 B k 2 
q 2 M k 2 + m 2 B ' 



As there is no transition, we expect that 

Tk 2 

G( k ) » ±f_ 

9m 



(15) 



(16) 



in the symmetric phase where is large. 

Deep in the broken phase, G(k) approaches zero, but 
causality implies that very long-wavelength (low k) cor- 
relations can only change slowly [THL IT^ | . We can give a 
rough upper bound for the rate of change, 



dlnG(fe) 



dt 

Using Eq. (fToT> . this becomes 

k 2 d In / 



<k. 



k 2 



dt 



<k. 



(18) 



Below the transition, lnmg w \/—m 2 /g 2 T, and if we 
keep on decreasing m 2 , then sooner or later Eq. JTSJl 
ceases to be satisfied for k less than some critical value 
k. The modes with higher k keep on decreasing and we 
approximate the final correlator by 



G(k) 



Tk 2 

exp 



2k 2 



(19) 



A Gaussian fall-off like this would follow naturally from 
diffusion, but our conclusions do not depend on the pre- 
cise form of the correlator, as long as it has a relatively 
sharp cutoff at k. The corresponding monopole density 
is given by 



n M 



d 3 k 
(2tt) s 



G{k) 



1/2 



(20) 



which agrees with Eq. Ijl3|) if we identify k = 1/£b. 

We can also find the monopole-monopole correlator 
in coordinate space by taking the Fourier transform of 
Eq. 0, 



G(r) 



Tk 5 e~ r k I 2 



it; (2 



and it is indeed positive at distances r< ^/Z/k. 

As a concrete example, let us now estimate the 
monopole density produced in the GUT phase transi- 
tion using only causality to limit the growth of It 
is clear that causality leads to a freeze-out, because the 
currect magnetic screening length would be proportional 
to exp(mjvf/2T) ~ exp(10 28 ) and therefore enormously 
longer than the size of the observable universe. This is 



still an oversimplification and the estimates should not 
be taken literally. 

At high temperatures, the effective mass parameter of 
the theory is m 2 (T) « g 2 (T 2 — Tq UT ) . Because of the ex- 
pansion of the universe, the temperature is decreasing at 
the rate dT/dt « -T 3 /M P , where M P « 10 19 GeV is the 
Planck mass. Near Tqut, we can therefore approximate 



2 £gut ^ 
y Mp 



(22) 



Deep enough in the broken phase, the monopole mass 
grows as 



m M 



t 9 

,2171 GUT- 



9 2 M P - 



(17) From Eq. we see that the growth rate of £b is 



dt 



T, 



11/4 



GUT 



gmJ^Mp 



lAf/2TGUT 



Tgut 
gM P 



x 7 ' A e* 



(23) 



(24) 



where we have introduced the dimensionless variable 
x = mM/2TcuT- We require that this is equal to 1 for 
the freeze-out scale, and find x ~ ln(gMp /Tqut), and 
consequently £b ~ 9 2 Mp/Tq VT . Then, Eq. ljT3|) tells us 
that the monopole density is 



n M 



1 



rpll 

J GUT 

Mp 



1/2 



(25) 



which we can compare with the prediction of the Kibble 
mechanism under the same circumstances |29j, 



Kibble 



g 2 T { 



GUT 



Mp 



(26) 



The two results differ by a factor of g 6 (Af P /T G uT) 3/2 , 
which is not particularly large for realistic GUTs, but 
could in principle have any value. 

According to Eq. (|12|) . the typical number of 
monopoles in a cluster is 



Af nct — 
IV M — 



Qm 



M P 



(27) 



This combination is, again, of order one, which means 
that there is a possibility of forming small clusters. 

As was already mentioned, the estimate in Eq. I|25|l is 
not very precise. The main factor in this is that the mag- 
netic charges are likely to move diffusively rather than at 
the speed of light. The true freeze-out scale £s is nec- 
essarily shorter than our estimate and therefore Eq. I|25() 
can be thought of as an approximate lower bound and 
Eq. I|27l) as an upper bound. Furthermore, if the tran- 
sition is fast enough, which may actually be the case in 
the GUT transition, the approximation in Eq. Q that 
the monopoles are point particles is not justified and one 
should instead use a field theory description. Neverthe- 
less, this simplified calculation shows the places where 
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more accurate physical input is needed to improve the 
estimates. 

It is also interesting to apply this same picture to 
cases where the GUT symmetry is restored only briefly 
after inflation, either because of "non-thermal" fluctu- 
ations 0, IE EJ or because the reheat temperature 
is slightly above Tqut- The estimated monopole den- 
sity depends on the low-momentum behaviour of G(k) 
given in Eq. ((TrJ|) . Because of charge conservation, the 
monopoles and antimonopoles must be produced in pairs, 
and even if they move at the speed of light, the leading 
term in G{k) grows as G(k) ~ nu^t 2 . It will therefore 
take at least the time t eq » (^T/njvr) 1 ' 2 w £s to achieve 
the form 1)16(1. This conclusion can also be reached by 
considering the time it takes for the pairs to reach the 
equilibrium size ~ £g. 



This means that if the GUT symmetry is restored only 
very briefly, for a period shorter than t Gq w (5 2 Igut) , 
the number density of monopoles will be suppressed. In 
reality, the equilibration process is probably significantly 
slower, and therefore t oq can be much larger, perhaps 
even so large that the bounds on the reheat temperature 
disappear completely. In any case, a more careful anal- 
ysis of the dynamics is needed to estimate how strong 
the suppression is in practice and whether it solves the 
monopole problem in the case of non-thermal symmetry 
restoration. 
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